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CsJ Abstract 
> ! 

Q I It is generally accepted that neutron stars form in core collapse events that are accompanied 
by a supernovae (types II or lb or Ic). Typical progenitors are, therefore, larger than ~ 2.1M0. 
We suggest ^1 12] that the binary pulsar J0737-3039 provides evidence for a new formation channel: 

^ I collapse of a light progenitor. This binary pulsar J0737-3039 has several remarkable features including 

^ I among others: a very tight orbit with a Keplerian velocity of 600km/sec, a low eccentricity, and a 

00 ' location 50pc from the Galactic plane implying that the system has, at high likelihood, a small 

^ ', (compared to Keplerian) center of mass velocity. A significant mass loss during the formation of the 

T-H I second pulsar would have lead either to an eccentric orbit or to a large center of mass velocity or to 

^ ] both. Therefore, we can set a strong upper limit on the progenitor's mass. A progenitor more massive 

^ I than 1.9Mq is ruled out (at 97% confidence). The kinematically favored option is of a progenitor mass 

JpH| around 1.45Mq. Recent evidence for a rather low velocity proper motion supports this prediction and 

O I decreases the likelihood of the standard (high mass progenitor) scenario. Lack of variations in the 

I pulses' profiles (which indicate no significant geodetic precession) provides further support for the no 

' kick and low progenitor mass formation scenario. 



^ ■ 1. Introduction 



The remarkable binary system J0737-3039 O EI was discovered during a pulsar search carried 
out using a multibeam receiver at the Parkes 64-m radio telescope in New South Whales. It is 
composed of two pulsars J0737-3039A and J0737-3039B denoted here A and B. This discovery provided 
a new ideal general relativistic laboratory. The eclipses of A beyond B provide a superb way to explore 
pulsar magnetospheres. The relatively short life time of this system to gravitational radiation emission 
has lead to a revision of the binary merger rate in the galaxy. We have suggested, immediately 
following the discovery P, that the orbital parameters of this system and its location close to the 
Galactic plane pose strong limits on the origin of this binary system and on the mass of the progenitor 
of the younger pulsar B. We review those arguments here and examine how do these predictions look 
now almost two years after the discovery. 

We begin by summarizing the relevant parameters of this system: The separation, R, (the 
sum of the semi-major axes) of the pulsars today is 8.8 • 10^'^cm. The orbit is almost circular and 
the eccentricity, e, is 0.087779. The periods Pa,b and their time derivatives provide upper limits for 
the life times of the pulsars: tA ~ 210Myr and ~ 50Myr. Dispersion indicates that the system is 
600pc from Earth. It is located 50pc from the Galactic plane. 

The present evolution of the system is determined by gravitational radiation emission. The 
separation of the pulsars and the eccentricity decrease with time and the system will merge approxi- 
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mately 85Myr from now. An interesting feature of the evolution is that the eccentricity decays with 
a faster rate than the orbital separation [S]. Thus an eccentric orbit becomes first circular and then it 
decays. We can also integrate backwards in time and determine the system's parameters at the time 
of its birth. We find that 50Myr ago, when the pulsar B was born, the separation was lO^^cm. The 
ellipticity at that time was only slightly larger e ~ 0.11 than the present one. Thus, the system was 
born with a low ellipticity. The values are not that different at 200Myr. At that time, the eccentricity 
was 0.14 and the separation was 1.2 • lO^^cm. As these values are not that different from the present 
ones our analysis is insensitive to the exact age of the pulsar. 

The binary B1913+16 was studied for almost 3 decades since its discovery jHj. Wex et al. 
[Zj used the ample data to place interesting limits on this system. For example, they found that 
the natal kick must have been directed almost perpendicular to the spin axis of the neutron star 
progenitor. On the other hand, the rather large post-SN eccentricity implies that the pre-coUapse 
orbital separation (I.8-4.6R0) or the mass of the progenitor (4-32M0) cannot be tightly constrained. 
The system J0737-3039 is different from previously detected binary pulsar systems in that it has a 
very low eccentricity, and it was small also after the collapse. Thus, the orbital separation before the 
collapse, for example, can be constrained to have been A = 1.0 ± 0.1 ■ lO^^cm = 1.45 ± O.ISR©. As 
the system J0737-3039 was much tighter than B1913+16, the orbital Keplerian velocities (and hence 
all other relevant velocities) are much larger. This would play an important role in our analysis. The 
location of J0737-3039 only 50pc from the Galactic plane suggests that the system has a small center 
of mass velocity. This is another important factor. 

The formation of the second pulsar is described according to current picture by a core collapse 
event that involves a supernova and mass ejection from the system. Dewi and van den Heuvel [Hj 
considered a scenario in which the progenitor star lost most of its envelope through interaction with 
its companion A. Prior to the formation of the second pulsar, tidal interaction between the progenitor 
and the neutron star has led to a circular orbit. The lost mass took place via a common envelope 
phase, at which point the companion J0737-3039A was spun up and its magnetic field was suppressed 
by accretion. They estimate that the progenitor mass was more massive than 2.3Mq. According 
to their estimates for lower masses the common envelope phase would evolve too fast and would 
result in too small separation. The formation of the second pulsar is described according to the 
current picture by a core collapse event that involves a supernova and mass ejection from the system. 
Standard evolutionary scenarios, lead neither to neutron star formation nor to core collapse, from 
progenitors that are less massive than 2.1 — 2.3M0 llOj. 

We begin with a general discussion of the orbital parameters of a binary after an instantaneous 
mass ejection event, such as the one that takes place in a supernova. We then apply these general 
arguments to J0737-3039. We show that that if the formation of the second (younger) pulsar involved 
significant mass ejection then the hole system would have had have a large center of mass motion, of 
the order of the large Keplerian velocity. However, with a large center of mass motion it is unlikely to 
find the system in the the Galactic plane. We use these arguments to put limits on the mass ejection 
and on the progenitor's mass. These arguments were put forwards right after the discovery of the 
system and they have led us to predict that the system will have a very small (tens of km/sec) proper 
motion velocity jl]. Later on upper limits on the velocity were found and we consider the implications 
of these upper limits to our considerations. We conclude with a discussion of the implications of these 
results to stellar evolutionary scenarios and to the rate of NS mergers. 
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2. Mass Ejection and Orbital Motion 

We consider first the influence of mass ejection during tfie formation of the second pulsar, on 
the orbital motion. We consider in the following a system that due to the tidal interaction between 
the older neutron star A and the progenitor of B was in circular motion prior to the second supernova. 
During this supernova a mass, Am, is ejected from the star B while the star itself becomes a neutron 
star. The thrust of the ejected mass. Am, gives a velocity, Vcm, to the center of mass (CM) of the 
remaining system. In addition, the mass loss would lead to an elliptic orbit or even to the disruption 
of the system. For a spherically symmetric mass loss, Vcm will be: 



where vk = \jG{mA + rriB)/ R is the Keplerian velocity of the two stars relative to each other, just 
after the explosion, with R being the distance between the two stars at that time. Within the context 
of J0737-3039, m^ = 1.377(5)M0, m^ = 1.25O(5)M0, Am is, of course, unknown and m^j = m^+Am 
is the initial mass of B while vk ~ 600km/s. With m^ ^ rriB, Vcm would be of the order of w/</2 
unless Am <^mB- 

The common explanation for a birth of a system with a low eccentricity (with Am m^) is for 
B to have had a natal kick jJJ . While the origin of these kicks is not clear, independent evidence for 
such kicks arise from the peculiar motion of individual pulsars ^2]- Surely, with suitable adjustment 
of the kick velocity the flnal orbit can have a low eccentricity and even become circular. However, in 
order to decrease the eccentricity the kick velocity given to B must be opposite to its actual velocity. 
This can be understood intuitively. After the mass ejection the mass of the system and hence the 
gravitational attraction between A and B decrease and the system flnds itself with "too much" kinetic 
energy. Hence the orbit is elliptic or even disrupted. To reduce the ellipticity we need to slow down 
B (there is no way to slow down A at this stage) and this can be done with a kick velocity opposite 
to its original velocity. 



Fig. 1. The kinematics of mass ejection and the orbital motion with and without a kick velocity. The 
conditions after a kick velocity that leads to a circular orbit are marked with dotted lines. 




(1) 





Even without a kick velocity the mass ejection causes the CM of the system to move in an 
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Fig. 2. The maximal kinematically allowed mass loss Animax and progenitor mass niBi for a given change 
in the CM velocity Avcm- The upper mass depends on the actual eccentricity (0.088 < e < 0.14), hence 
the finite width of the line. The hatched area is kinematically excluded. The region allowed by standard 
evolutionary scenario through the formation of a He-star, which requires a minimum mass of 2.IM0 
(marked A ^01) to 2.3M0 (marked B [H]) to form a NS through a core collapse SN, is marked as shaded 
triangles. The no-kick solution that results with the above range of eccentricities is marked by the short 
heavy line. 



increases Vcm further and the above value can be considered as a lower limit. In other words, the 
nearly circular orbit today implies either a large CM velocity, roughly as given by Eq. ^ or a small 
ejected mass Am ^ m^. 

Fig. 2 depicts change in the CM velocity (relative to the unknown CM velocity prior to the 
explosion) for a system with two masses moving on a circular orbit, assuming the new binary system 
attains the initial eccentricity of J0737-3039 (0.088 < e < 0.14). In the following we assume that 
this initial CM velocity was small and we approximate Vcm ~ ^t^cm- We will return to this point in 
the conclusions. Fig. 1 shows that the minimal CM velocity for a 2.IM0 progenitor is 120km/s. The 
no-kick solution marked here is different from the one of Dewi and van den Heuvel [HI who assume 
that A rather than B has undergone the second collapse. 

3. Limits on the CM Motion 

The position of the system almost in the Galactic plane sets an additional constraint (to the 
one imposedby the low eccentricity) and it has led us to predict shortly after the discovery of 
the pulsar that it would have a very low CM velocity and hence a low peculiar velocity . Since that 
time there have been several different estimates of the CM velocity of this system. Ransom et al. 
[T3] have estimated Vcm±, the CM velocity of the binary on the plane of the sky, using the observed 
scintillations of the system. They find a rather large value: Vcm± = 141 ± 8.5km/s (with 96.0 ± 3.7 
km/s along the orbit and 103.1 ± 7.7km/s perpendicular to it). This value excludes the region in Fig. 
1 left of a vertical line of ~141km/s. These findings were questioned recently by Coles et al. ^1] who 
suggested that the scintillation pattern is anisotropic. When including anisotropy, they find a much 
lower value: VcmX = 66 ± 15km/s. Pulsar timing give even lower upper limit on the CM velocity of 
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Wcm± < 30km/s The last two values are actually consistent with each other as the estimate of 
66km/s was not corrected for the motion of the Earth (TSj. The region in Fig. 1 to the right of the 
vertical line of ~30km/s is consistent with this observation. 

The observed distance of the system from the Galactic plane, enables us to place a statistical 
upper limit on Vcm- Stars move in a periodic motion in the vertical direction. For small vertical 
oscillations, the potential of the Galaxy is harmonic: $ = 27iGpoz'^, where po ~ O.25M0/pc^ is the 
mass density in the disk ^H]- This gives a vertical orbital period, Pz ~ 50Myr. The typical velocity 
for an object at Zo^s is ~ 27iZobs/Pz- ^obs ~ 50pc implies then that the expectation value of the 
vertical velocity is of the order of 6km/s. 

To quantify the probability for having a particular CM velocity given the observed Zobs we 
perform Monte Carlo simulations that follow the formation of the system. We assume that star B 
had a given mass rriBi, and that a randomly oriented kick Vkick was given to it. We also assume 
that the progenitor distribution has an initial Gaussian distribution in the amplitude of the vertical 
oscillation, with a width = 50pc (other cr2<100pc gave very similar results). At the moment of 
formation, we assume it had a random phase within its vertical motion. We calculate the CM kick 
velocity v^m, and assign it a random direction, then integrate the vertical motion of the pulsar for 
50Myr using a realistic Galactic potential |17j . 
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Fig. 3. The probability that a binary system will end up within 50pc from the Galactic plane, with 
0.088 < e < 0.14 and with a transverse velocity less than 30 km/s, given that SOMyrs before, the 
progenitor system had a circular orbit, that star B had a progenitor mass ms^i, and it obtained a random 
kick velocity of size t'kick- The probabilities are normalized to the most likely Msi, vi^^ck- The He-star 
solution requires rriBi larger than about 2.1M0. With a fine tuned Vkick) this type of a solution is ruled out 
at the 99.7% c.l. Solutions with niBi ~ 1.5ib0.2MQ and Wkick^SOkm/s are kinematically more favorable. 

Fig. 3 depicts the probability that a system with a given niBi and Vkick could find itself with 
0.087 < e < 0.14 after 50Myr, within 50pc of the Galactic plane and with a transverse velocity less 
than 30km/s as given by pulsar timing ^Hj- Other probability distributions with different assumptions 
on the current status of the system can be seen in O . With the constraint that the transverse velocity 
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is less than 30 km/s, a fine tuned rriBi = 2.IM0 model can be ruled out at even 99.7% c.l. Even 
without this constraint on the velocity, without any information on the CM velocity, we find [2] that 
even if V]^ick is fine tuned to be near either 130 or 315km/s, a system with a mass of msi = 2.1Mq 
could result with the observed configuration in only about 3% of the random realizations (as compared 
with the most favorable conditions having lower m^i and Wkick)- Other t^kk's, or higher mass systems 
are kinematically even less likely. On the other hand, low ejected mass solutions are favored. This 
was the basic argument of p. 

An inspection of Fig. 3 reveals that there is a kinematically favorable solution with a small 
mass loss and natal kicks ranging from to ~ 50km/s. Since a small mass loss necessarily implies 
a small natal kick, the solutions with Vkick<^50km/s are physically unlikely. If we therefore limit 
ourselves to Wfcjcfc^30km/s, then without fine tuning, a large fraction of the progenitor systems will 
result with the observed configuration with a progenitor mass of 1.45MQ<mBj<1.65M0. Thus, a 
mass loss of about 0.3 ± O.IM0 is most probable. The results are qualitatively the same even if we 
do not enforce that the resulting system has a CM velocity of 66km/s, or if we replace the condition 
on the eccentricity to e < 0.14 rather than 0.088 < e < 0.14. 

Willems and Kalogera fHl estimated the allowed Am and Vkick assuming, following standard 
scenarios, 2.1 < rriBi < 4.7M0. They find 50 < Vkick < 1560km/s and a most likely value of 150km/s. 
However, in their analysis Willems and Kalogera implicitly assume that the current radial CM velocity 
of the system (which has not been measured) is large, impliying that the system is moving with a very 
large velocity almost directly towards us. In fact their assumption on the magnitude of this velocity 
dictates their resulting probability distribution of the kick velocities. Our calculations include the 
geometric probability that a large CM motion will be directed towards us. When this factor is included 
we find, as explained above that the "standard" high mass progenitor and large kick scenario is ruled 
out. 

4. Conclusions and Implications 

We are left with two physically distinguished scenarios for the formation of pulsar B. In the 
first, the progenitor is a kinematically "unlikely" but theoretically plausible 2.1 — 2.3M0 He-star 
progenitor — around the minimal masses estimated from stellar evolution scenario |Hl 13 ^] • In the 
second scenario, the progenitor is a kinematically favorable ~ I.5M0 young stellar core. The more 
probable low mass solution requires a new type of stellar collapse. Intermediate solutions are neither 
statistically favored nor do they fit any plausible theoretical scenario. 

Our conclusions were based on the low eccentricity of the observed system and on its location 
near the Galactic plane. When we first suggested this model in January 2004 the CM velocity of 
the system was not known and we have predicted a low CM velocity p]. This was confirmed later 
by pulsar timing observations ^H] that have shown that the system has a very low peculiar motion. 
An independent supportive evidence for out low (or no) kick model is given by the fact that pulsar 
A did not show variations in its pulse shape in fifteen months of observations fH]- Such variations 
were expected [20] due to geodetic precession. While this observations may results from a particular 
precession phase, it is more likely that it arises due to alignment of the spin of pulsar A and the orbital 
angular momentum. This suggests, in turn, that the orbital angular momentum hasn't changed in 
the second collapse and indicates a low kick velocity during the formation of B*. 

* A large kick velocity would presumably have a component out of the original orbital plane and would have changed 
the direction of the orbital angular momentum. 
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Before turning to the implications of these two solutions we consider, first, three assumptions 
made in our analysis, (i) We have assumed that prior to the formation of the second pulsar the 
system was in a circular motion. This follows from all evolutionary scenarios that lead to a neutron 
star and a small mass progenitor (even a 2.3M0 is a very small mass progenitor), (ii) We have 
assumed that the second mass loss was instantaneous, namely, shorter than a fraction of an orbital 
period. Given that the orbital motion is of several hundred km/s while typical mass ejection velocities 
in SNe are higher than 10,000km/s, this assumption is reasonable (for all conventional neutron star 
formation scenarios), (iii) We have assumed that the CM velocity, prior to the formation of the second 
pulsar was small. One would expect that the system would have acquired a CM velocity during the 
formation of the first pulsar. About half of the system's mass was lost during this event and this 
should have resulted in some CM velocity. However, this velocity would have been of the order of 
half the Keplerian velocity of the system at that time (see Eq. Q). Given the fact that the orbital 
separation was much larger we could reasonably expect, but not prove, that this velocity was of order 
of several tens of km/s. Moreover, if the system would have acquired a large CM velocity in the first 
SN, there would have been an even smaller probability to find it in the Galactic plane today. 

We turn now to the most likely scenario, the very low mass scenario. We imagine the same 
evolutionary scenario in which some time before 50 Myr, system A and progenitor B were in a 
common envelope phase and B lost most of its mass keeping practically just its core of ~ 1.45Mq. 
This progenitor leads to a small CM motion and does not require any kick velocity. This solution is 
kinematical preferred. However, it requires a new mechanism for the formation of the pulsar as He- 
stars of I.45M0 do not collapse to form neutron stars jHlQD]. The observation that the progenitor mass 
is very close to the Chandrasekhar mass leads us to conjecture that the process involves the collapse 
of a supercritical white dwarf. For example, the progenitor may have been a degenerate bare core just 
above the critical Chandrasekhar mass, supported by the extra thermal pressure against collapse. As 
it cooled, the additional support was lost and the core collapsed to form a neutron star. A second 
possibility is that it was formed just below the Chandrasekhar mass, and as it cooled, neutronization 
at the core increased the baryon to electron ratio, and with it reduced the Chandrasekhar mass until 
the progenitor became unstable and collapsed. Note that the object must have been composed of O- 
Ne-Mg as a collapsing CO core would have carbon-detonated and it would have exploded completely 
forming a type I SNe and leaving no remnant. However, it seems that the parameters space for these 
two scenarios is quite small. A second possible formation scenario is that the progenitor of B had a 
core of about 1.4M0 and it was gradually stripped of its envelope until it eventually collapsed as a 
bare core with almost no mass ejection. In either case it is clear that whatever the formation scenario 
was it is drastically different from the "standard" scenarios for NS formation. 

This new solution passes an immediate non trivial test. With a small mass loss, only small kick 
velocities are possible and in fact we can estimate in this case (see Fig. 3) the mass loss needed to obtain 
the initial eccentricity e ~ 0.11. We find (while conservatively assuming that the collapse took place 
SOlso^Myr ago) Am = e{l + q)mB = O.28±O.O7M0 and a progenitor mass of 1.53±O.O7M0, which is, 
indeed, just above the Chandrasekhar limit. Intriguingly, some mass loss, in the form of u losses of a 
few times lO^^ergs, must take place. The estimated Am corresponds to E^, ^ Amc^ ^ 4.2 • lO^^ergs, 
which is in the right range. Of course, if some mass is ejected as well, E^, will be smaller and a small 
kick could arise, but the total mass-energy lost will be the same. The consistency of this mass and 
energy loss with the previous physical picture increases our belief in this new and unusual scenario. 

It is intriguing that (after the publication of our first suggestion P) other arguments has led 
Van den Heuvel 121] to conclude that about half of the pulsars that are in binary systems had low 
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mass progenitors. 
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